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Abstract 

It is well known that for low speed flows the use of the compressible 
fluid dynamic equations is inefficient. The use of an explicit scheme 

requires At to be bounded by 1/c. However, the physical parameters change 
over time scales of order 1/u which is much larger. Hence, it is not 
appropriate to use explicit schemes for very subsonic flows. Implicit schemes 
are hard to vectorize and frequently do not converge quickly for very subsonic 
flows. We shall demonstrate that if one is only interested in the steady 
state then a minor change to an existing code can greatly increase the 
efficiency of an explicit method. Even when using an implicit method the 
proposed changes increase the efficiency of the scheme. We shall first 
consider the Euler equations for low speed flows and then incompressible 
flows. We then indicate how to generalize the method to include viscous 
effects. We also show how to accelerate supersonic flow by essentially 
decoupling the equations. 
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Euler Equations for Subsonic Flow 

We first consider low speed flows for rotational inviscid flow. Since 
the flow may be rotational we consider the Euler equations rather than the 
potential equation. We only consider schemes in conservation form. The use 
of conservation form allows the same code to be used for highly subsonic, 
transonic and supersonic flows. 

The Euler equations, in two space dimensions, can be expressed as 


w^ + f + 
t x 


g y = ° 


( 1 ) 


where (x,y) represent general curvilinear coordinates. Since we are only 
interested in the steady state we replace (1) by the system 


M * w. + f + g = 0. (2) 

t x y 

The requirements on M are that the matrix be nonsingular and that the 
resultant equations be well-posed. It is straightforward to solve (2) with an 
explicit scheme. With an implicit method only the diagonal portion of the 
matrix to be inverted is changed. Though the code solves (2) we shall only 
analyze the constant coefficient problem 

M * w„ + Aw + Bw = 0 (3) 

t x y 

where the matrices M, A, B are constant. Let w^) = Tw, Aq = T A T 

Bq = T B T” 1 , M” 1 = T M”* 1 t” 1 , where T is chosen appropriately [6], then 

(3) can be converted to 


I 


- 2 - 




q = Y u - X v, r = X v - Y u (6) 

n y y ’ xx 

q and r are the contravariant components of the velocity and (X,Y) are 
the Cartesian coordinates. 

We now consider the case of low speed flows. We wish to choose so 

that the eigenvalues of Mq Aq and Mq Bq are independent of c. We also 
wish to choose Mq to be positive definite. This will imply that (4) is a 
symmetric hyperbolic system and so well-posed. One choice is 

c 2 

— 0 0 0 

z 2 

0 10 0 

0 0 10 

0 0 0 1 

where z 4- = max(e , u 4 ’ + v 4- ). e is introduced so that the matrix M is not 
singular at stagnation points. Transforming back to curvilinear coordinates 
we find that M = I + dQ, d = (y-l)(z 2 /c 2 -l)/c 2 
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where 

s 2 = (u 2 + v 2 )/2, h = c 2 /(y-l) + s 2 . 


( 8 ) 


(9) 


We note that given the first row of Q the following rows are derived by 
multiplying the first row by u, v, h respectively. Hence the product of 
Q times a vector requires only six multiplications. 

o o o 

Let M = z /c then the largest eigenvalue of D = A sin 0 + B sin <(> 
is given by 


1 _ |w|(l + M 2 ) + X/ 2 (l - M 2 ) + 4(a 2 + b 2 )z 2 
A = ^-2— 


( 10 ) 


where 

w = q sin 6 + r sin <(> , a = Y sin 0 - Y sin <{>, b = X sin <f> — X sin 0. 

y x x y 

We see that near a stagnation point M = 0(e), X = 0(e). While at M = 1, 

X = | w | + /a 2 + b 2 c. Hence, at small Mach numbers the largest eigenvalue, 
and hence the time step, is independent of the sound speed c. At transonic 
sound speeds the largest eigenvalue is comparable to the nonconditioned 
case. Hence, the preconditioned problem allows a larger time step for all 
subsonic flows (see also [5] for the isoenergetic case). 
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We next consider 
given by 


incompressible flow. In conservation form the 


u + v = 0 


x y 

u + (u 2 + p) + (uv) = 0 
u x y 

v t + (uv) x + (v 2 + p) y = 0 


system is 


(id 


where p is the pressure normalized by the density. We wish to integrate 
this system using the artificial compressibility method [1]. If p t /c 2 is 

inserted in the first equation the system is hyperbolic but not well 
conditioned [2]. Instead we replace (11) by 


— j p + u + ( u 2 + p) + (uv) = 0 
c x y 

P t + V t + (uv) x + (v2 + p) y = °> 


( 12 ) 


c is a given nonzero function. It is evident that at the steady state both 
systems coincide. The new system can readily be shown to be unitarily 
equivalent to a symmetric hyperbolic system and so is well-posed and well- 
conditioned. The extra time derivatives can be considered as a precondition- 
ing matrix similar to that previously considered. It remains to decide how to 
choose the function c. For c larger the coefficient of many time deriva- 
tives is small, however the time step of an explicit scheme will be large. 
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Hence, in this case c is merely a scaling of the time scale. Instead, we 
wish to choose c as large as possible without the need for decreasing in any 
substantial way the stability criterion. In one space dimension the stability 
criterion for (12) is given by 

f~ 2 2 

At < u + /u + 4c (13) 

Ax 2 


Hence, a reasonable choice for c is 


2 

c 


max 




+ V 



(14) 


Viscous Flow 

We next consider the incompressible steady state Navier-Stokes equation 

u + v = 0 
x y 

(u^ + p) + (uv) = 4* ( u + u ) (15) 

v y R ^ XX yy' 

(uv) + (v^ + p) = 4‘( v +v )• 
v 'x v v/ y R ^ xx yy J 

As before we shall consider a pseudo time dependent approach to the steady 
state* We consider two ways of extending the previous results to include 
viscous effects* 

The first possibility is to use a leapfrog method for the inviscid part 
and a Dufort-Frankel scheme for the viscous portion [1]. The inviscid part 
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72 2 

bounds the time step by Ax//u + v while the viscous part is 
unconditionally stable. The second possibility is to use a semi -implicit 
method based on the preconditioned system. The advection terms are treated 
explicitly while the viscous terms are treated implicitly. As before, the 
explicit part restricts the time steps inversely proportional to the 
velocity. The implicit part only contains linear terms. Hence, the implicit 
part can be inverted without any need for an iterative method. 


Supersonic Flew 

We next consider supersonic inviscid flow. In this case the matrices A 
and B can be simultaneously diagonalized by a congruence transform [4]. 
Hence we can find a matrix M so that (2) is similar to a diagonal system 
(for the linearized constant coefficient problem). Thus, In the supersonic 
regime we can choose a different time step for each equation. The resultant 
system is still hyperbolic and well-posed. 
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